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CL|' Abstract 

D : 

^ • A number of DR renormalization constants in softly broken SUSY-QCD are eval- 

uated to three-loop level: the wave function renormalization constants for quarks, 
^SJ ■ squarks, gluons, gluinos, ghosts, and e-scalars, and the renormalization constants for 

^ I the quark and gluino mass as well as for all cubic vertices. The latter allow us to 

■ derive the corresponding (3 functions through three loops, all of which we find to be 
identical to the expression for the gauge /3 function obtained by Jack, Jones, and 

■ North [1] (see also Ref. [2]). This explicitly demonstrates the consistency of DRED 
I with SUSY and gauge invariance, an important pre-requisite for precision calculations 

■ in supersymmetric theories. 

a^ ■ 
o 



1 Introduction 

Up to now, particle collider experiments could not provide clear evidence for physics 
beyond the Standard Model (BSM). The quest for extensions of the Standard Model comes 
either from non- accelerator observations, such as Dark Matter and Dark Energy, or it is 
based on purely theoretical considerations. For example, naturalness [3l[31[5] requires some 
mechanism to stabilize the Higgs mass at the electro- weak scale (fine tuning problem); the 
intriguing evolution of the gauge couplings towards a common value at high energies needs 
an enlarged particle spectrum if unification is to occur at scales compatible with proton 
decay O [71 [8] ; any explanation for electro- weak symmetry breaking presumably requires 
an embedding of the Standard Model into a higher symmetry. 

In the light of these considerations, supersymmetry (SUSY) is a strong candidate for an 
extension of the Standard Model [9]. Not only does it provide solutions for Dark Matter, 
it also solves the fine tuning problem by canceling the quadratic divergences in the Higgs 
self energies, it explains electro-weak symmetry breaking by a simple evolution of the 
parameters of the Higgs potential, and unification of the gauge couplings can be achieved 
for rather natural choices of the parameters. 
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This latter issue is a rather unique feature: it actuaUy allows one to make quantitative 
predictions at energy scales that are several orders of magnitude larger than what can be 
achieved at current and probably also at any future particle collider. Crucial ingredients 
for such indirect measurements are precision data as well as precision calculations (see, 
e.g., Ref. [ID]). 

In order to study the unification of the gauge couplings in SUSY, one needs to include the 
effect of the super symmetric particle spectrum into the running. In addition, dimensional 
regularization with minimal subtraction, called MS, is no longer a good renormalization 
scheme in the sense that it explicitly breaks supersymmetry. An alternative regular- 
ization, so-called Dimensional Reduction (DRED) was suggested by Siegel [11], however, 
inconsistencies of this method were pointed out only shortly afterwards [12]. Neverthe- 
less, renormalization by combining DRED with minimal subtraction (the DR scheme) has 
become the preferred scheme in higher order supersymmetric calculations [13]. A thor- 
ough formulation of DRED which isolates the source of the inconsistencies has been given 
in Ref. [Tl]. It was pointed out that, although a mathematical consistent formulation of 
DRED is possible, it could violate SUSY at higher orders of perturbation theory. It is the 
main goal of the current paper to demonstrate this consistency within SUSY-QCD up to 
the three-loop level. 

Renormalization group functions, governing the energy dependence of masses and cou- 
plings, are among the simplest quantities to compute in perturbative quantum field the- 
ory. For example, the anomalous dimensions of the strong coupling constant in QCD 
(the QCD P function) and the quark masses are known at four-loop level in the MS 
scheme \15\ \TE[ \T7\ [T8] (the conversion to the DR scheme was done in Ref. [19] ) . 

For supersymmetric gauge theories, one can devise a particular renormalization scheme, 
so-called NSVZ [20], where an all-order relation between the gauge (3 function and the 
anomalous dimension of the chiral supermultiplet is valid. So, in the absence of the matter 
supermultiplet, i.e. for a SUSY Yang- Mills theory, the /3 function is known to all orders 
in the coupling constant. The problem with this scheme is that general conversion rules 
to schemes used in perturbative calculations such as MS or DR are not known. However, 
by explicit calculation of the abelian /3 function in the DR scheme, comparison with the 
NSVZ result, and arguments based on its holomorphy, the authors of Ref. [1] found the 
two-loop conversion formula for the coupling also in the non-abelian case. The SUSY-QCD 
(3 function in the DR scheme is thus known to three loops. The result was later confirmed 
by an explicit calculation using the background field method [2] , which requires three-loop 
corrections to the gluon propagator. 

Applying the same method based on holomorphy of the NSVZ scheme to softly broken 
SUSY gauge theory, the authors of Ref. [21] derived the renormalization group equation 
governing the running of the gaugino and squark masses, valid to all orders in the coupling 
constant. Using this result, an elegant relation between the gaugino and the gauge /3 
functions was formulated in Ref. [22], both for the NSVZ and the DR scheme. 

One of the goals of the current paper is to add another confirmation of these results by 
calculating for the first time directly the relevant vertex diagrams through three-loop order. 
Dimensional reduction is implemented by explicitly introducing e-scalars, the (4 — d)- 
dimensional components of the gluon field. 
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In addition, we derive the gauge (3 function I3s from all possible three-point functions in 
SUSY-QCD that involve the strong coupling a<j- The fact that in each case we obtain 
the same expression provides an explicit check on the consistency of DRED with gauge 
invariance and super symmetry. In particular, we calculate the renormalization of the 
quark-squark-gluino vertex which is not fixed by gauge invariance but by supersymmetry 
alone. Furthermore, we confirm the relation between (3s and the gluino mass anomalous 
dimension 7^ by an explicit calculation of 7^ to three loops. We also verify the three-loop 
formula for the quark anomalous dimension 7^ |23] that was derived in the superfield 
formalism. It is an essential ingredient for the derivation of the four-loop SUSY-QCD 
gauge P function in the approach of Ref. [23] . 

The e-scalars occurring in DRED have Feynman rules that can be derived from the ones 
of the gluon by replacing d-dimensional objects with 2e-dimensional ones. In a supersym- 
metric theory, the coupling constants for the e-scalar and the gluon must be identical, 
which means that the /? function for the e-scalar couplings must be given by the one of the 
strong coupling constant a^. We check also this relation through three loops, by evaluating 
three-point functions involving e-scalars through three loops. 

The remainder of this paper is organized as follows: In Sect. [21 we review the method of 
DRED and describe its implementation in our setup. Sect. [3| describes some details of the 
calculation for the Green's functions and the (3 functions, and also gives the main results. 
Concluding remarks are given in Sect.[H 



2 Dimensional Reduction 
2.1 Notation and technical setup 

We implement DRED by introducing the quasi-four, -d, and -2e-dimensional spaces [24] 
(see also Ref. [11]) Q4, Qd, Q2e, where d = 4 — 2e and 

Q4 = Qd'JQ2e and QdriQ2, = 0. (1) 

Also, we introduce quasi-/our, -d and -2e-dimensional sets of indices, implying the corre- 
spondence 

zy^(z>,p), p^{p,p), a^{a,a), (2) 

in the sense that, for example, 

t^' = t^^ + etc. (3) 

for any vector t. I.e., and denote the projections of G Q4, to the subspaces Qd and 
Q2e, respectively. For generic indices without reference to (quasi-) dimensionality we will 
use the letters a, /3, 7, 6. 

For the metric tensor, we have 

... . - (4) 

<7'^% = ^^ 5'^% = ^, g^^ = 4, gl = d, gl = 2e. 
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The usual commutation relations for four-dimensional Dirac matrices are carried forward 
to their d- and 2e-dimensional projections: 

{7^7'^} = 2/^ {7^7''} = 2/^ {^(^,y} = o. (5) 

These relations together with the constraint Tr 1 = 4 allow us to evaluate fermion traces 
(without 75) in the following way: 

Tr(7^i . . . ^Pi . . . ^Pm) ^ ^Tr(7'^i • • • 7'^")Tr(7'^i • • • 7^^'") . (6) 

However, the quark-squark-gluino vertices involve 75, and the problems involving 75 in 
other than four dimensions are of course well known. We treat it as follows: 



1. Use the relations 



{7^75} = {7^75} = o, (75)' = !, (7) 

until there is at most one 75 per fermion trace. If all 75 matrices are eliminated in 
this way, one can continue with Eq. 

2. Fermion traces with a single 75 matrix are reduced to the form Tr(7"7^7'^7''75) by 
using Eqs. ([5]) and ([7]), and then evaluated by applying the formal replacement 

Tr(7"7^7^7''75) = 4z e"^T^ . (8) 



3. The intrinsically {not quasi-)four-dimensional Levi-Civita symbol 

with quasi-dimensional indices makes sense when combined with a second one of its 
kind. By applying appropriate projectors, our calculations are set up in such a way 
that this is always the case, so we can use 

^afS-yS^ _ Jo' „P „7 „<5] /q\ 

e Ea'p'-^'S' — 9[a' 9 9y 9 s'] > 1^) 

where the square brackets denote anti-symmetrization. Note that due to g^P" = 0, 
one can write, e.g., 

e"^ ef,'u>p>a' - 9ijx> 9 0>9 p>] 9a> , ef,>o'p'a' - 9^^, 9 o'\ 9[pi 9 , 

e^^'^^eA'^'s'r' = , etc. 

We made sure that the terms that need to be treated in this way generate no higher than 
simple poles in e. Therefore, we can be sure that the above procedure is valid without 
introducing additional finite counterterms. Of course, the fact that we find the same /3 
function from all vertices also supports our treatment of 75. 
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2.2 e-scalars and corresponding Feynman rules 



We work in = 1 supersymmetric QCD with n/ quark flavors, corresponding to the strong 
sector of the MSSM. The Feynman rules are therefore given by the usual ones for QCD, 
involving quarks q, gluons g, and ghosts c, plus the Feynman rules for their supersymmetric 
partners: the squarks qi, qn, and the gluinos, g. In our implementation, instead of the 
partners of the left- and right-handed quarks, we will use the mass eigenstates 



The mixing angle 9q depends on the squark masses; however, all the quantities in this 
paper will be independent of any masses, and therefore, they must also be independent of 
the squark mixing angle. This serves as a welcome check on our calculations. 

The tree-level Feynman rules for SUSY-QCD in DREG can be found in the literature; in this 
paper, we will use those given in Ref. [25]. DRED results are obtained by supplementing 
these Feynman rules by new ones for the so-called e-scalars, as will be explained in what 
follows. 

DRED is defined by assuming that all fields depend only on the coordinates of a 
(i-dimensional subspace of the usual four-dimensional Minkowski space. I.e., the 2e- 
dimensional components x^' have no physical relevance. In particular, this means that 
derivatives and momenta can always be replaced by 9^ and p/x. The technical im- 
plementation is done by decomposing the quasi-four-dimensional vector fields into d- and 
2e-dimensional components: 



For convenience, will be simply called the "gluon field" in what follows, while will 
be explicitly referred to as "four-dimensional gluon field". Afi is the so-called e-scalar. 
The corresponding Feynman rules for the e-scalars are obtained from the d-dimensional 
SUSY-QCD ones by replacing the d-dimensional indices by 2e-dimensional ones. In this 
way, one arrives at the following vertices: qqe, gee, gge, eeqq, ggee, eeee, where e denotes 
the e-scalar. Their Feynman rules are given explicitly in Appendix lAl 

The bare coupling constants for the gluon and the e-scalar are identical by construction. 
However, a priori it is not clear whether this holds also for the renormalized couplings. 
In principle, all e-scalar couplings could be different without violating gauge invariance. 
In fact, it is well known that they differ in standard (i.e., non-SUSY) QCD (see, e.g., 
Refs. |26 1 1191 [27] ). Even more: in order to renormalize the quartic e-scalar vertex, one has 
to take into account all possible color structures for it, and attribute to each one a separate 
coupling constant. For SUSY-QCD, we have explicitly checked that at the one-loop order 
only the (3 function associated with the usual color structure of the four-gluon interaction, 
i-e. fabefcd^, docs not Vanish and it equals the one-loop gauge /3 function. So, through 
one-loop, one can identify the coupling constant of the corresponding e-scalar quartic 
interaction with the strong coupling constant and set to zero the other three couplings. 

^fabe denotes the structure constants of the gauge group. 




Afj, — Ap_ -\- A^ . 



(10) 



5 



(a) (b) 

Figure 1: (a) Feynman diagram involving Majorana fermions in a non-trivial way. Solid 
lines are quarks, dashed lines are squarks, slashed springy lines are gluinos, and the exter- 
nal lines are gluons. The arrows on the lines denote the charge flow. Note that the fermion 
line involves quarks with opposite charge flow, (b) The same diagram after applying the 
prescription proposed in Ref. [28^ I29j. The arrows on the gluino lines denote the external 
flow (see text); lines with two "clashing" arrows are "flipped" quarks which have external 
flow opposite to the charge flow. 



This order of accuracy is sufficient for the results discussed in this paper, as the e-scalar 
quartic interactions contribute to the anomalous dimensions starting from the two-loop 
order. 

In order for the renormalized Lagrangian to obey SUSY, the decomposition of Eq. (jlOp 
should hold also at higher orders. Therefore, the renormalized gluon and e-scalar coupling 
constants must be equal, i.e., their P functions must be the same. We will explicitly 
verify this for the qqe, the gge, and the gee vertex through three loops, thus confirming 
consistency of DRED with SUSY at next-to-next-to-next-to-leading order of perturbation 
theory. 



2.3 Majorana fermions 

Another technical issue is given by the Majorana character of the gluinos. A well-defined 
and practical prescription how to deal with Majoranas in higher order perturbative calcu- 
lations was given in Ref. [28\ [29] . Unfortunately, using this method in combination with 
the Feynman diagram generator qgraf [30] is not straightforward. We implement it by 
parsing the output of qgraf with a PERL program, giving each chain of fermion lines a 
well-defined direction. If the chain contains a Dirac fermion / whose ordinary fermion 
direction ("charge flow") is opposite to the one chosen by the PERL program ("external 
flow"), we replace this particle with a "flipped fermion" / which has opposite charge flow. 
Majoranas have no definite charge flow and are simply interpreted as Dirac particles with 
charge flow equal to the external flow. 
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llg 

2lg 
3,1 
1,4 
5,2 
2,6 
4,3 
5,3 
6,4 
6,5 



Sql , sql 
fgm,fgm 
fgm,fgm 
Sql , sql 



fQ,fq 

fQ,fq 
fQ,fq 
fQ,fq 



llg 
2lg 
3,1 
1,4 
2,5 
6,2 
4,3 
5,3 
4,6 
6,5 



fQ,fq 
fQ,fq 
f Qx,fqx 
f Qx,fqx 
Sql , sql 
fgmx,fgin 
fgmx,fgin 
Sql , sql 



Figure 2: The diagrams of Fig. [T] in topological notation. 



An example is shown in Fig.[TJ In topological notation, diagram (a) is represented by 
the left column of Fig. [21 where each line denotes a propagator in obvious notation. Af- 
ter application of the prescription of Ref. [28l [29], the diagram in Fig.[T](b) results, or, 
again in topological notation, the right column of Fig. [3 One can see that the quarks 
(fQ,fq) connecting vertices 2 with 5 and 6 with 2 have been replaced by "flipped" quarks 
(fQx,fqx), and the Majorana notation for the gluinos (fgm,fgm) has been changed to the 
Dirac notation, where particle and anti-particle are distinguishable (fgmx.fgm). 

The result looks like regular qgraf output with only Dirac fermions. The flipped fermions 
have Feynman rules that are easily derived from the original fermions using Ref. |28l [29] . 
They are given explicitly in Appendix [Al Finally, the overall sign is determined in the 
usual way by counting the number of closed fermion loops. 

Note that the same approach has already been applied successfully to earlier calcula- 
tions [3T] . 

For the remainder of the calculation, we use the following setup: the diagrams are brought 
to MINCER notation [32] with the help of the C++ programs q2e and exp [33]. The actual 
calculation is done in FORM |34j : After applying appropriate projectors in Lorentz, Dirac 
and color space and evaluating the fermion traces, we use MINCER to evaluate the Feynman 
integrals in terms of a Laurent series in e. 

The representation of the four-gluon vertex through an auxiliary, non-propagating particle 
allows us to factorize the color factor from the rest of the diagram and to evaluate it 
separately. For that purpose, we use the program color. frm [35] which expresses the 
result in terms of color invariants. 

3 Evaluation of the renormalization constants 

We closely follow the method of Refs. [361 [3Z]) we calculate connected, amputated 
Green's functions and derive the coupling constant renormalization from Slavnov- Taylor 
identities. For example, if we compute the A^-point Green's function with external fields 
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(^1, • • • ,(t)n and denote its coupling constant by g, we obtain 

5- /7 (11) 

where the Z^. are the wave function renormalization constants for the (pi, Z^^...^^^ is the 
corresponding vertex renormahzation constant, and Zg the charge renormahzation. 

Since we work in a minimal subtraction scheme, the renormalization constants (and thus 
also the (3 function) are independent of any mass scale. We remind the reader that due 
to this fact, decoupling of heavy particles has to be done "by hand" in the DR scheme, by 
introduction of decoupling constants for the renormalized parameters of the model. The 
decoupling constants in the MSSM for the strong coupling constant and the quark mass 
have been calculated to two- loop order in Refs. [TOt [38]. 



The anomalous dimension for the coupling g can be obtained from the simple pole of Zg 
in the usual way via 

where g^ is the bare coupling, and /i is the (arbitrary) renormalization scale. Zg does not 
explicitly depend on /i, and if g is the only coupling constant in the theory, then 

d —eg 

Applying DRED to a non-SUSY theory like standard QCD, Zg in general depends also on 
the e-scalar couplings, and Eq. (jl3p assumes a slightly more complicated form (see, e.g., 
Ref. [26]). The P function is defined through 

P{as) = fi -T-^ , where = — . (14) 

d^^ TT 47r 

In analogy to Eq. ()14p we define the anomalous dimension for fermion (quark and gluino) 
masses 

^2 drui din Zm,i . , /, 

7i = — = -'^Ps—. , i = q,g, (15) 

nii d/i^ das 

where Zm,i is the mass renormalization constant. 

Since we work in a SUSY theory, we expect the renormalization constants for the e-scalar 
couplings to be identical to the ones of the strong coupling ^f^, provided DRED is consistent 
with SUSY. When evaluating the strong coupling renormalization Zg at n-loop, at most the 
(n — l)-loop expression for the e-scalar coupling renormalization Z^ is needed. Therefore, 
the calculation of Zs at n-loop order checks the identity Zg = Zf. through (n — l)-loop 
order: if Zg 7^ Z^. at (n — l)-loop order, in most cases the n-loop expression for Zg would 
not be local. However, in some cases it can also happen that simply the 1/e pole is wrong. 

In our calculation it is n = 3, so we check the equality to two-loop order. In addition to 
that, however, we have explicitly evaluated a number of e-scalar couplings to three-loop 
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7^ loops 


1 


2 


3 




7^ loops 





1 


2 


3 


cc 


1 


14 


423 




ccg 


1 


2 


77 


3920 


qq 


4 


86 


3583 




qqe 


1 


6 


319 


21669 


ee 


7 


100 


3902 




qqg 


1 


8 


439 


30078 


99 


6 


130 


5577 




99^ 


1 


8 


445 


31815 


qq 


8 


157 


6760 




qqg 


1 


6 


430 


36868 


99 


11 


171 


6954 






1 


14 


618 


38741 










999 


1 


12 


657 


46314 










qqg 


1 


12 


674 


52205 










ggg 


1 


26 


1105 


70705 



Table 1: The number of diagrams contributing to the Green's functions evaluated in this 
work. Left table: two-point functions; right table: three-point functions. The first column 
indicates the external legs of the Green's function, the other columns show the number of 
diagrams at the individual loop orders. The number depends on the implementation of 
the quartic vertices. In our calculation, we split the gggg, ggee, and the eeee each into two 
cubic vertices by introducing non-propagating auxiliary particles. 

order (see Table [T]) and also find equality to the strong coupling constant. In this way, 
consistency of DRED with SUSY is confirmed through three- loop order. 

For the calculation of the renormalization constants within the minimal subtraction scheme, 
one is free to choose any masses and external momenta, as long as infra-red divergences 
are avoided. We set all masses to zero, as well as one of the two independent external 
momenta in the three-point functions. We therefore arrive at three-loop integrals with 
one non- vanishing external momentum q which can be calculated with the help of MINCER. 
As a check, we have also calculated some of the three-point functions with non-zero par- 
ticle masses m by expanding them in the limit m? /(f <C 1 with the help of asymptotic 
expansions [39] • In the final expression the limit m ^ can be taken and the result coin- 
cides with the one obtain with the massless set-up. Possible infra-red singularities would 
manifest in \nm? /q^ terms which are absent in our calculations. 

Traces with a single 75 do not contribute to any of the two-point functions that we calcu- 
lated. They do contribute for some of the three-point functions though, in particular the 
qqg, the gge, and the qqe vertex. An example diagram for the latter vertex is shown in 
Fig. [31 Such diagrams contribute (among others) a color factor d'^'^'^d"^'^'^ (for the notation, 
see Refs. [351 127]). but they cancel against the same factors from other sources in the final 
result for the renormalization constants and the (3 function. 

The results for the renormalization constants can be obtained in electronic form from 
Ref. [iQ]. It is straightforward to obtain from them the three- loop SUSY-QCD (3 function 
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Figure 3: Sample diagram for the three-loop qqe vertex where a non-vanishing trace with 
a single 75 matrix occurs. The internal line styles are as in Fig.[TJ 



which reads 



(3i 



3 



Ca- 
3 ^2 



21 



E 

n>0 
1 

2 



V vr 



(16) 



1 



13 



16 



16 



16 



where Cp = (?^c ~ 1)/(2't-c)i = ?^c are the Casimir operators of SU(nc), and 2Tj = nj 
is the number of quark flavors (which is equal to the number of squark flavors). The result 
in Eq. (|16p is in full agreement with Ref. [U [2] . 

The case Tf = (SUSY Yang-Mills theory) has been treated before, both in the super- 
field formalism [23] and the diagrammatic approach [19]. Full agreement has been found 
between the two methods up to four-loop order. 

The anomalous dimension for the quark mass is independent of any SUSY-breaking pa- 
rameters and is given to three-loop level by 



7g(a,) = - ^ 



n>0 



n+1 



In 1 



lt = -\c'F + lcACF-\TfCF, 



3 

16 



1 



1 



- + -C(3) + - - -C(3) CaCf 



16 



(17) 
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where C is Riemann's zeta function with ^(3) = 1.20206. . .. The result of Eq. (jl7p is in 
agreement with Ref. [23]. 

In Ref. [22] it has been shown that the anomalous dimension of the gluino mass defined 
through 

7s(«.) = -5^(^) ^ In: (18) 

n>0 

can be related to the gauge (5 function of Eq. (fTUj) via 

7^ = (n + l)(3n ■ (19) 

We have explicitly verified that this relation holds through three loops for the case of 
vanishing SUSY-breaking parameters except the gluino mass. 

4 Conclusions 

The main result of this paper is a confirmation of the consistency of DRED with SUSY at 
three- loop order. 

We have shown that the same result for the SUSY-QCD (3 function is obtained from all 
three-particle vertices, even the ones involving e-scalars which are used for the implemen- 
tation of DRED. Furthermore, we verify a predicted relation between the /? function and 
the gluino mass anomalous dimension to three-loop order and confirm the result for the 
quark mass anomalous dimension present in the literature. 

Along with this paper we provide the results for all two- and three-point renormalization 
constants [ID] which should be useful for other calculations within the MSSM. 

Acknowledgments. This work is supported by DFG through SFB/TR 9 and contract 
HA 2990/3-1, and by the Helmholtz Alliance "Physics at the Terascale". We would like 
to thank K.G. Chetyrkin, D.R.T. Jones, P. Kant, and D. Stockinger for illuminating 
discussions and useful comments. 

A Feynman Rules 

In this appendix we give the Feynman rules involving e-scalars and flipped fermions. The 
notation is the same as in Appendix A. 2 of Ref. [25], except that all Lorentz indices oc- 
curing there are understood as d-dimensional ones. Note that in that paper, the sign of 
the ggqq vertex should be reversed, and there is a typo for one of the qqg vertices which 
should read 
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r « 
a 




igsT-,V2{nlPL-n1^PR) 



r 



a 




igsT,%V2{nlPR-n'^,PL) 



Also, the four-gluon vertex did not contribute in Ref. [25j . so for completeness, let us quote 



it here: 
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5 ^ 


^ 9 



■ 2 f abe ^cde 
-Wsf J 



au ha up i^cr 

9^ 9^ - 9^^ 9 



The following are the tri- linear vertices involving e-scalars (all momenta incoming): 
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As introduced in the main text, q denotes a "flipped" quark whose "external flow" (indi- 
cated by the arrow on the hne) is opposite to its charge flow. For clarity, the external flow 
has also been indicated on the gluino lines in the egg vertex. The other Feynman rules 
involving flipped quarks are given as follows: 





igsT-,V2{nl^PL-nlPR) 



r 

q ► 



a 




9 



igsT-,V2{n1,PR-n1,PL) 



The quartic vertices involving e-scalars are 



9 

f), b 
9 



,- e 

~~ e 



ade fbce 



jace joae _j_ jaae j 



/Z, a \ / c, -p -ig^J'^'^r"' [g~^Pg'~" - g'^^'g 
p ft da - ialr^f"' {9^'9'~" - 9^~"9' 
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